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NOTICE 


This report was prepared as an account of Government- 
sponsored work. Neither the United States, nor the National Aero¬ 
nautics and Space Administration (NASA), nor any person acting 
on behalf of NASA: 

a. Makes warranty or representation, expressed or implied, 
with respect to the accuracy, completeness, or usefulness 
of the information contained in this report, or that the use 
of any information, apparatus, method, or process dis¬ 
closed in this report may not infringe privately-owned 
rights; or 

b. Assumes any liabilities with respect to the use of, or 
for damages resulting from the use of any information, 
apparatus, method, or process disclosed in this report. 

As used above, "person acting on behalf of NASA" includes 
any employee or contractor of NASA, or employee of such contractor, 
to the extent that such employees or contractor of NASA, or employee 
of such contractor prepares, disseminates, or provides access to, 
any information pursuant to his employment with such contractor. 

Requests for copies of this report should be referred to: 

National Aeronautics and Space Administration 
Office of Scientific and Technical Information 
Washington 25, D. C. 

Attention: AFSS-A 
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1. 0 INTRODUCTION 


This report contains all the results obtained to date with the use 
of a computer simulation of the tape recording process, the tape 
demagnetization-remagnetization cycle, the readout process and several 
detector processes. 

After a description of the programs (JPL-4 and JPL-6), the report 
is divided into two parts. The first deals with the pure tape recording studies 
made. These include all phenomena up to and including the reproduce head 
output. 

The second part deals with the problems associated with detection, 
drop-outs, noise and other signal processing operations. 
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PRECEDING PAGE 3LANK NOT FILMED. 

2. 0 DESCRIPTION OF THE FOURIER SERIES COMPUTER 

PROGRAM (JPL-4) 


2. 1 Program Deck JPL-4 
2. 1, I General Description 

This deck simulates a tape recorder in the recording and reproduction 
of digital information, and the detection of the reproduced signal. The 
program is designed to operate with a periodic input signal in which each period 
contains a six-bit binary digital information in a self-clocking code. The 
program computes first the magnetization recorded on tape (neglecting magnetic 
particle interaction) then the voltage output of the reproduce head, and finally 
the detector output. 

2. 1. 2 Input-Output 

Table 2.1 gives a full description of the data cards needed with this 
program. A set of these five data cards is needed for each reproduce head 
output voltage curve that is computed. If only such an output voltage curve 
is needed, these five cards, with FLAG set to 0.0 is sufficient. If, however, 
the detector output must be computed, then a zero sample signal and a one 
sample signal must be computed, as described below. Each of these sample 
signals requires in turn the computation of its own output voltage curve. 

When detector output computations are made, the first five data cards, with 
FLAG set to 1.0, provide the zero sample signal, and the next five data cards, 
with FLAG set to 2.0 provide the one sample signal. Following the compu¬ 
tation and storage of these two sample signals, each additional set of 5 data 
cards, with FLAG set to zero will result in the computation of a reproduce 
head output voltage curve, the cross correlation functions of this curve with 
respect to the zero sample and the one sample, and the detector output. 
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Table 2. 1 Data Card Description for JPL-4B Program 


Data 

Card No. 

Format 

j 

Variable 

Type 

Description 

1 

8F10. 0 

RLAMBDA 

REAL 

The length along the magnetic tape, in microinches, of one 6-bit period of the input signal. 



G 

REAL 

One-half of the record head gap width, in microinches. 



HSAT 

REAL 

Saturating value of magnetic field, in oersteds. ’ 



HMIN 

REAL 

Switching value of magnetic field, in oersteds. 



. HO 

REAL 

Deep-in-gap magnetic field in the record head, in oersteds. 



A 

REAL 

Head-to-tape spacing, in microinches. 



DR 


Sum of the head-to-tape spacing and tape thickness in microinches. 



SYM 


If the number of ones in the 6-bit code is even, SYM =1.0; if this number is odd, 

SYM = -1.0. 

2 

8F10. 0 

FLAG 

REAL 

If the reproduced signal for a 6-bit code is to be computed, FLAG =0. If a zero sample 
is to be computed, FLAG = 1.0. If a one sample Is to be computed FLAG = 2.0. 

3 

; 8110 

NPP 

INTEGER 

The number of points used to represent the magnetization on each half-bit of tape. 



NH 

INTEGER 

The number of harmonics used to represent the magnetization on tape. 



NY 

INTEGER 

The number of values of y (see Fig. 1) at which the magnetization is computed. 

I 

4 and 5 

8F10, 0 

PULSE 

REAL 

PULSE is an array of 12 numbers (either ones or minus ones) that represents the 6-bit code. 
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The program now provides for printing out the input data and the 
detector output. The user can easily provide for additional print-outs as 
needed. 


In different applications of this program it is desirable from time 
to time to plot a wide variety of different computed variables. One might, 
for example, wish to plot the magnetization on tape, the reproduce head 
voltage either or both of the cross ccMrelation functions, or the detectability. 
The detailed operations needed for plotting are most easily provided by a 
special subroutine. Plots can then be provided at any desired points in the 
program simply by inserting at these points the appropriate call statements 
to the plot program. 

2. 1. 3 Computation of Magnetization on Tape and Output Voltage from 
Reproduce Head 

Figure 1 shows the dimensions used in the formula given below. The 
distance above the record head to'the point at which the magnetization is 
computed is Y and the distance along the tape axis from the record head center 
line to this point is X. This is the axial distance to that point from the point 
on tape that is above the record head center line at t = 0. Thus 

X = X+vt (1) 

r 

where v is tape velocity and t is time. The component of record head magnetic 
field in the X direction is H . 

X 

The magnetization on tape is first computed using a non-interacting 
model. Demagnetizing effects are considered later, as indicated below. The 
remanent magnetization on the tape is assumed to be entirely in the longitudinal 
direction. The magnetization M, written at the point (X, Y) on tape, is given 
as a function by H by 

X 
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Coating v 
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\a ('Wx Cx,'^. f)) 


(2) 


where t is the point in time at which the point (X, Y) passes through 
o 

the second head center plane and t^, 1 ^ i ^ n are the points in time 
between successive bits. In the last term on the right of Eq, (2), n is 
the total number of non-zero values of M [ H , (X, Y, t.)] and X is 
expressed in terms of X^ and t^ by Eq. (1). In turn, 




o 



if 

1 Ha 

if 

V 


(3) 




(«- -"■) 

(4) 

if 

W 1 


(5) 

and 

■ 

K 

1 


(6) 

if 

1 1 

> 



where x is the susceptibility of the magnetic material (see Fig. 2)^ 
is the threshold value of the magnetic field at which switching can take 
place and is the saturating value of the magnetic field. In turn, 
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Fig, 2 Non-Interacting Model of Tape Remanence 
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(3 t~lx 

TO 


(7) 


( 8 ) 


(9) 


( 10 ) 


The record head field longitudinal component is 


H. rx,’y,f)= 


Ho 

11“ 




( 11 ) 


where H (t) is the deep gap field and g is the gap width, 
o 
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The deep-gap field, H^(t) is a periodic function of time, with 

a period of twelve self-clocking bits. For example, H (t) for twelve bits 

o 

of the pattern 0, 0, 1, 1, 1, 0 is shown in Fig. 3. Note that if b is the time 
interval of one bit, then for a pattern that contains an odd number of ones (for 
example that shown in Fig. 3) 




( 12 ) 


and for a pattern that contains an even number of ones 


N„(t) 


(13) 


Next the alogarithm simultaneously takes into account the effect 
of magnetic particle interaction, or demagnetization and computes the 
reproduce head flux. The first step toward this objective is to compute an 
average magnetization M given by 

di 




aid 


d‘ 




(14) 


where a is the head-to-tape spacing and d is the effective thickness of the 

magnetic material on tape. * Following this, we represent M by 

a 


no 


M 


a 


(Xy) « ao + 

o *1 


*For discussion of effective thickness of tape, see section 3. 5, 
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6 bits 


6 bits 



Fig. 3 Input Current Versus Time 
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where 



air 


This Fourier series representation is possible because M is periodic in 

cl 

12 by Eq. (2) through Eq. (14), In addition, from Eq. (12) and (13), 
and certain other previous equations. 




if the 6 bit code contains an odd number of ones and 


if the 6-bit code contains an even number of ones. As a result, M contains 

’ a 

only even harmonics if the 6-bit code contains an even number of ones 
and M contains only odd harmonics if the 6-bit code contains an odd 
number of ones. 


Mallinson^ has shown that the flux in the reproduce head $(t) 


is given by 


0rt)= 00 










(16) 


, 


where 



(17) 



(18) 


^ J, C. Mallinson, "Demagnetization Theory for Longitudinal Recording, " 
IEEE Trans, on Magnetics, Vol. MAG-2, No. 3, September 1966. 
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In Eqs, (17) and (18) 

[)^= fy)kc/') ^ COS )“i (19) 

j , rvt fc^ 0 ^ cqtK ('vile,/') -tj/'.r I >^k fviI< c/ % <• d fy\ M1 ^ /^q) 

M,AA 2 cosk rv\kc|')+ fv\k(3/')) 

y(lkyip>^U(nk^ ')(^a ^ (kK'>V4 

^ st^s (^k 

and ■ VI 

where g^ is the reproduce head gap width. 


The output voltage from the reproduce head is given in terms of the flux by 


V/ft) = W 


cjd)(t) 

c/ t 


( 20 ) 


where N is the number of turns on the reproduce head. From Eqs. (1), (16) 


and (20) 


VJO 


m _ 

y\^ * 



(21) 


The reproduce head output voltage is computed by Eq. (21). 
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2. 1. 4 Computation of Sample Signals and Cross-Correlation Functions 

To compute a zero sample, we first find the output voltage from the 
reproduce head for a signal input to the record head that consists of six zero 
bits. Since for this pattern 


Then 

and 


H (t + B) = -H (t) 
o o 

K (tH) = -V. (f) 


( 22 ) 


From Eq. (22), we see that there must be a time t at which V is zero, and 
’ or’ 

that 



From Eqs. (23) and (24) 

TJt) = S„(ft.OU = 0 

and from this and Eq. (25), we see that S is continuous at i at t + 2b, 

’ o (J o ’ 

as shown in Fig. 4 . 
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Fig. 4 Sample Signals 
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To compute a one sample, we first find the output voltage from 
the reproduce head for a signal input to the record head that consists of 
six one bits. Since for this pattern 


the 


and 


/-/„ (t* 

^4) = V^.(t) 


(26) 


From Eq. (26) there must be a time t^ at which is zero, and that 




(27) 


where m is any integer. To construct the one sample,, S^(t) we let 

1, s-Z&t-l-l, 

(28) 

t< t. t- 

(29) 

From Eqs, (27) and (28) 

S', (t,) - S, Ct,4 - o 

and from this and Eq. (29), we see that is continuous at t ^ f-ndt^ t-3/2b. 

The cross-correlation function of the signal with resp''''Ct to the 
zero sample, (t) and with respect to the one sample (i) aco given 
in this application by 
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C Cr) 

C,(y)- 

where X(t) is the signal. In the computer program, the integrals in Eqs. (30) 
and (31) are computed using the trapezoidal rule. 

Note that t^ and t^ are each chosen from a multiplicity of values. 

Furthermore, from Eqs. (30) and (31) the phase relationship between 

C (t) and C, (t) changes if 
o 1 


is changed as a result of a new selection of t^^ or t^ or both. The program 

has a provision to keep this phase relationship from changing, and to 

insure that S (t) andg, (t) are orthogonal. After S (t) and S,(t) are computed by 
o 1 o i 

Eqs. (26) through (29), S^(t) is slid along the time axis so that its first at 
absolute maximum occurs at the same value of t as the first absolute maximum 
of S (t). Then S (t) and S,(t) have the time relationship shown in Fig. 3. 

O O 1 

We define the "detectability", 6 (t) to be 

fCY). - I CoCr]l ,32) 




X(t)JAt*T) it 

"t ' t* U? Is 

xct)y,('t-r) i-t 


( 30 ) 


(31) 
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and compute 6(t) from t = 0 to t = 6b. Then the detector performance is 
computed using the following algorithm: 


1. Let T = 0 
o 


2. For i = 1, 2y.6 find such that 

|,SrY;)/ ^ I §■ Cr)l 

< r < r./H 


(33) 


3. Let S + be the set of integers i for which 

JTfrj >c>< 

and S - be the set of integers i for which 

S'cn) 


/ ^ I 

where Os " ^ 




r 




J 


t^O 


Xs ^ k) 




J. 




(34) 


Then compute 

P and Q from 


P = Vvicvy 
• 


(35) 

A 1 h 


0 - ;< 

1 ^ (rJI 

^^6) 

) 1 K 



The quantities 

P and Q are then used in error rate ci Uri '! 
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Another parameter required for the error rate compilations is p, defined 
as 

and computed by this equation. Finally, the ratio P/p and Q/p are 
computed and printed out 

2. 1. 5 Comparison of Variable Names in Text and Computer Program 

Table 2.2 shows the variable names in the computer program 
opposite the corresponding symbols from the text. 

2. 2 Description of the Computer Error-Rate Calculation (JPL-6) 

2. 2. 1 General Description 

This deck is used to compute error rates in digital detection. The 
error rate is computed on the basis of the detected zero and one signal 
P and Q, (computed in JPL-4 B), the noise power per unit bandwidth, N^, 
and the parameter jS, defined by 

fQ-JO) 

pCa) --P ^ 

where a is the head-to-tape spacing, and p(a) is the probability that the 
head-to-tape spacing will have the value "a". 
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Table 2. 2 Comparison of Varible Names in Text and in 
Computer Program JPL-4B 


Name of Program 

Name in Computer 

Symbol in Text 

Number of Equation in 

or Subroutine 

Program 


which Symbol is Used 

Main Program 

A 

a 

14 


ALPII 

a 

34 


ALPHA 

n 

15 


C 

a 

n 

15 


CROSS 

C^(r) 

3 0 


C ROSSI 

c ^(t) 

31 


D 

a + d* 

14 


DFA 

D 

n 

19 


G 

g/2 

1 1 


HO 

H 

o 

1 1 


HMIN 


3 


RESX i 

! 

M (x ) 
a r 

14 


RFSXZ 

v^(t) 

21 


RESX3 

1 

6{r) 

3Z 


RHO 1 

P 

3 7 


RK 

k 

lb 


B 

b 

n 

I 5 


SAMPLE 0 

S (t) 
o 

24 


SAMPLE 1 

S^(t) 

28 


X 

X 

r 

21 


Y 

y 

11 

Subroutine 

H 

H 

1 1 

MAGCAL3 

HO 

X 

H 

o 

1 i 


RES 

M(x^, y) 

2 


X 

X 

11 


Y 

y 

1 1 

Function 

A 

a 

19 

DEMAG 

D 

a + d' 

19 


DEMAG 

D G 
n n 

19 


GR 

g/2 



RK 

k 

19 


RMUl 

", 

19 


RMU2 

"2 

19 

Subroutine 

CROSS 

6 (t) 

32 

DETECT 

PPLUS 

P 

35 


PMINUS 

Q 

35 


PP 

P/p 



PM 

Q/p 



RHO 

P 

37 
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2. 2. 2 Input - Output 

Table 2.3 shows the format of the data cards used for this program. 

The deep-in-gap magnetic field and the number of bits per inch (both which 
are input) are not used in the computation but appear in the printed output as an 
aid to identification of the computed error rate values that are printed out. 

For each error rate computation, the computer prints out, on a single line, the 
values of N^, H^, the bits per inch, and finally, the error rate. The 
program prints out the values of N^, /3, A (head-to-tape spacing), P/yo, and 
Q/p, in the order in which they are read in. To permit the computation of a 
number of error rates, the program returns to read data card No. 4 (in Table 
2.3) at the completion of each computation. 


2. 2. 3 Error Function Complement 


The algorithm for error rate computation employ the computation 
of the error function complement. The error function complement is computed 
in the subroutine named ERROR. This subroutine uses algorithms, provided by 
Abramowitz and Stegun. ^ If the argument exceeds 2. 5, the algorithm of 
paragraph 7. 1. 26 of this reference is used; otherwise, the algorithm of 
paragraph 7. 1. 23 is used. 


2. 2. 4 Computation of Error Rate 

The quantity ER(a) is computed by 



(39) 


where the ratios P/p and Q/p were computed by means of the program deck 
JPL-4B. Note that ER(a) is a function of a because P and Q are functions of a. 
Finally, the error rate, ER is computed by 


^Milton Abramowitz and Irene A. Stegun "Handbook of Mathematical Functions, " 
National Bureau of Standards, Applied Mathematics, Series No. 55, June 1964 
pp. 298-299. 


RR 67 -36 


21 


AMPEX 



a-- Q, 

pCrO ER cJ a +-5 C 




CK^;io 


(40) 


where p(a) is given by Eq. (38), ER(a) is given by Eq, (39) and a^ is the 
head-to-tape spacing beyond which errors occur even in the absence of 
noise. The first term on the right of Eq. (40) is the probability that 
noise will cause error, and the second term on the right of this equation is 
the probability that error will result only from excessive head-to-tape 
spacing. This integration is carried out in the computer program by use of 
the trapezoidal rule. 


2. 2. 5 Comparison of Variable Names in Text and in Computer Program 
(JPL-6) 

Table 2.4 gives the names used in computer program JPL-6 for the 
variables used in the above discussion of the error rate computation. 

Table 2. 3 Data Card Description for JPL-6 
Computer Program 


Data 

Card No. 

1 

2 

3 

4 

5 

6 

7 

8 


Format 

8110 

8F10. 0 
8F10. 0 
8F10.0 

8110 
8F10. 0 
8F10. 0 

8F10. 0 


Variable 

NN 

NEX 

RNO 

EX 

HO 

BPI 

NA 

A 

P 

Q 


Type 

INTEGER 

INTEGER 

REAL 

REAL 

REAL 

REAL 

INTEGER 

REAL 

REAL 


Description 

This is the number of values of to use in the computation. 

This is the number of values if ^ (sec Eq. 38)) to be used in the computation. 
This is an array that contains the values of N^ to be used. 

This is an array that contains the values of ^ to be used. 

This is the value of the deep-in-gap field. 

This is the number of bits per inch to which the computation pertains. 

NA is the number of value head-to-tape spacing to be used in this computation. 
This is an array that contains the head-to-tape spacing used in this computation. 

P is an array that contains the values of P/p (see Eq. ( 39)) for the values of heac 
to-tapo spacing just read in. 

Q is an array that contains the value of Q/p for the values of head-to-tape 
spacing read in. 
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Table 2. 4 Comparison of Variable Names in Text and 
in Computer Program JPL6 


Name in Computer Program 

Symbol in Text 

A 

i a 

ERFCP 

erfc {P/{p V2N^ )) 

ERFCQ 

erfc (Q/{p )) 

EX 


P 

P/p 

PA 

P(a) 

PROD 

P(a) ER(a) 

Q 

Q/p 

RNO 

N 


o 

SUM 

ER 
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3. 0 PRE-DETECTION STUDIES 


3. 1 Absolute Output Voltage Multiplier 

With the sole exception of program JPL-2, the original Fourier 
transform deck, all the computer outputs are now subject to the same 
calibration factor on voltage. 

If the number of turns on the reproduce head is N and the maximum 
tape remanence is 4ffM^ gauss and the track width is (omils and the tape 
^^locity is V ips, then the correct output voltage is derived by multiplying the 
computer output by the factor; 

6. 25*10"^^ (N) (4ffM ) (to) (V) 

r 

For example if; N = 250 turns, 4 irM^= 1000 gauss to= 80 mils 
and V = 60 ips (these typical figures being close to IBM compatible practice), 

^^^true " ^^^calc 25* lO'^^ . 250 . 1000 . 80 . 60 


or 


(c) = (c) , X 7. 

^true calc 


-3 

10 volts 


Since (c)^^^^^ is approximately 1 unit, we have (c)^^ ~ 7. 5 milli-volts 
for this example. 
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As another example, let us consider an audio recorder of the type 

used in the calculation of the noise power figures* Here we will typically 

have; N= 1000 turns, 4irM = 1250 gauss, co = 50 mils, V = 10 ips, which 

^ -4 

yields a conversion factor of 4. 10 . The tape background noise level of 

-11 2 -11 "4 2 2 

10 units /cycle thus becomes 10 ( 4 .10 ) volts /cycle which 

2 

equals 160 (nanovolts) /cycle. 


Finally, it must be remembered that none of the programs include 
any allowance for the well-known flux "shunting" effects in reproduce 
heads. Since typical head efficiencies (i. e., the ratio of flux threading 
the output coil to total tape flux) are approximately 20 to 40%, the computer 
outputs may be anticipated to be proportionately high. 


3. 2 Comparison of the Fourier Series Program Output Voltage with 
Experiment 


A test was made of the ability of the six bit program (Fourier 
Series) to duplicate pulse crowding or bit interaction effects correctly. 
Typical results at 3200 BPI are shown in Fig. 5, where the reproduce head 
output voltage for the six bit NRZ work 111000 is plotted against tape 
longitudinal distance for three values of the deep gap record field 
(H^ = 1500, 2250 and 3000 oe). Other parameters in the calculations, 
which correspond closely with the actual case shown dotted, are: 


head to tape spacing 
tape coating thickness 
record gap length 
reproduce gap length 
range of switching fields 
permeabilities 


a = 20 p" 

d = 400 (x" 

g = 500 jj." 
rec 

g , = 250 li¬ 
re prod ^ 

H^/H^ = 100/500 oe 
— 5/2 


*See section 4. 
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It will be observed that all the essential features of the actual 
signal are reproduced by the computer. A small discrepancy is apparent 
in that the 2. 0 X SAT computer curve matches the actual curve better than 
the 1. 5 X SAT curve which should show the best match. The computer is 
in fact slightly underestimating pulse crowding. At BPI less than about 
2000, pulse crowding is due simply to the superpositicn of adjacent pulses. 

At BPI greater than 2000 as is the case here, additional pulse crowding, due to 
the interference of the record zones, occurs. Such crowding is , of course, 
non-linear and non-super-imposable, and is dependent critically upon the 
record gap length, record deep gap field and the specific record head field 
function being used. The present computation error is most probably due 
to the use of the Karlquist head field expressions. 

No additional work upon the problem will be performed during 
the present contract. It may be safely concluded that no worse match than 
that shown will occur, for any bit pattern, providing the BPI < 3000 for the 
400 |x" thick media or BPI < 6000 for the 200 |ji" media. 

A more exhaustive program of experimental validation was under¬ 
taken as an extension to the work reported here. The details may be found 
in Appendix II. Suffice it to say that, upon such close examination, the 
program was found to be useable to higher bit densities than indicated above. 

In the saturation record level case, the program is good to at least 5, 000 
BPI; at lower record levels, (which would be used) good agreements are 
found to be at least 10, 000 BPI. 
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3. 3 The Effect of Record Current Finite Rise Times 

A test was made of the effect of non-zero record current (or more 
pertinently record field) rise times. This investigation was made for 
the reason that the actual experimental waveforms, taken at a tape speed 
of 75 IPS, used a waveform generator with 0, 5 p sec rise time. The rise 
distance along the tape is thus 37 to which may be added perhaps 
another 30 jji" due to the limited frequency response (eddy currents) of the 
record head structure. The total rise distance of the record field is then 
probably about 7 0 p", which is an appreciable fraction of the bit length 
at 3200 BPI (315 p"). It was felt to be worthwhile to establish that the 
rise time effect was not contributing to the discrepancies observed. 

The Fourier Series program was modified to permit exponential 
behaviour of deep gap field (see Fig. 6 below). 



Fig. 6 Exponential Behavior of Deep Gap Field 
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The factor is simply the exponential time constant multiplied 

by the tape velocity. Because mathematical complications arise if the rise 

distance is long compared to the bit interval, the program is restricted to 

X ^ 1/2 X BPI. Figure 7 shows two waveforms produced for the six bit 
o 

word (NRZ) - 111000 at a deep gap field of 1500 oe with two rise 
distances, 0 and 156 p". The latter distance is one-half the bit interval* 
at 3200 BPI. 

It will be observed that, apart from a phase shift, wherein the 
finite rise distance case is moved upstream on the tape with respect to the 
zero rise distance case, there is indeed very little effect. It is concluded, 
then, that providing the rise distances are not sufficiently long to cause 
severe "staircasing" (i. e., incomplete switching) of the record field, the 
effects are negligible. 

This conclusion is not surprising, in view of our previous finding 
that the isolated pulse output waveform was virtually independent of the 
precise (saturation) recording conditions. The written transition relaxes 
(i. e., demagnetizes) and retains little information at wavelengths less than 
the coating thickness. A corollary of this conclusion is that, of course, 
nothing useful can be gained by any of the more or less elaborate schemes 
of pre-shaping the input record current. 
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3,4 Output Anomalies 

In the course of a series of over 100 runs investigating the peak- 
to-peak output voltage of an "all ones" pattern on a function of deep gap field 
(H^), BPI and head-to-tape spacing, it was found, whenever the record 
head-to-tape spacing was comparable with the bit interval being used, and 
whenever the record-head field was below the level required to saturate the tape 
(8. e,, 500 oe contour beyond the the back side of the tape), that the output 
voltage was an oscillatory function of the deep-gap field, (See Fig, 8 .) 

This oscillatory behaviour was not foreseen; neither could any 
reference to it be discovered in the literature. It was assumed, erroneously 
that some mistake in the computer program was causing the effect, and 
nearly ten days (and almost 20 computer runs) were spent examining a variety of 
possibilities, A slight computer program error was discovered (concerning the 
algebraic sign of the increments in remanent magnetization acquired by tape 
elements) but rectification of this fault did not stop the oscillations. It was 
finally found, by hand claculation, that the oscillations occur as a direct 
consequence of the noninteracting remanence loop model [ see MTPR #1, 
feq. 1)]. The phenomenon only occurs when several field reversals occur 
during the time (or distance) that a tape element is passing through the 
critical range of field strengths (H^ and H2) on the downstream side of the 
record head. 

It is not known, at this time, whether such oscillation really occurs 
in practical systems. A point worth noting is that it is predicted to occur 
only under conditions where the output signal is already at least 20 db 
(factor of ten) below the peak level possible. It may, therefore, have 
easily escaped observation. The whole phenomenon is closely related to 
the old hi-fi buff's technique of d. c. erasing tape with a "horse shoe" 
permanent magnet which is placed obliquely against the tape. In this case, 
the tape experiences a large longitudinal field, followed by a smaller 
field in the opposite direction, and can thus be left with zero remanence. 
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It is of interest that if the non-interacting M^-H model were to be 
replaced by an empirical model of remanence acquisition, the effect would 
still persist. 

It is concluded, then, that the effect calculated most probably 
would be found in practice but is not of any great consequence. It was 
decided, therefore, to leave the computer program unaltered. 

3. 5 Modification of the Program to Correct the Depth of Recording 

It will be recalled that after each lamina of tape has been recorded 
upon, the next step in the computer program is the computation of the average 
magnetization; 

_ , a + d’ 

M = — r M dy (See section 2. 1) 

d *’a 


The points of concern here are the value of the upper limit of integra¬ 
tion (a+d') and the effective tape thickness d' to be used when the whole 
depth of the coating is not completely magnetized (i. e., saturated). 

Previously, this effective tape depth (d') has been calculated according 
to the smaller of 

d* = d 

or a + d' = depth of 100 oe contour on head center plane. 

Thus any material which did not experience even the minimum switching 
field (100 oe) was merely ignored throughout the remainder of the calculation. 
This procedure obviously gives the maximum reasonable value for d*. The 
minimum reasonable value would be given by ignoring all material which is 
not saturated, i. e, : 
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d' = d 

or, a + d' = depth of 500 oe contour 

It will be appreciated that much of the remaining program is 
critically dependent (at low values of the deep-gap field H^) on the 
particular criterion chosen. This occurs mainly because the demagnetization 
and remagnetization processes are very sensitive to the product k d ( or k d') 
(see page 5, MTPR #3). The old criterion (100 oe contour depth) led to the 
pulse width being almost independent of deep-gap field (see Table IV, 
page 5, MTPR #3), because even at H = 500 oe, the 100 oe contour still 

7 7 O 

penetrates a 200 |jl" coating completely. The result is the over-estimation 
of pulse crowding at low record currents. 

Obviously no single correct criterion can be decided upon. This 
unfortunate situation is forced upon us by our desire to use the closed 
form solutions of the demag-remag problem. The alternative procedure using 
iteration is believed to be of virtually astronomical magnitude (200 x 
200 calculations per bit cell per iteration). 

It has been decided to resolve this problem by using an "engineering" 
approximation. The effective depth shall be considered to be half way 
between the maximum and minimum values possible; i. e., we take d' to be 
the smaller of 

d' = d 

or a + d' = depth of the 300 oe contour 

It is noted, in passing, that Kostyshyn (IEEE Trans. Mag, 2^ 3, page 
236, Sept. 1966) assumed the minimum reasonable depth. In other words, he 
ignores completely the partially magnetized material. 
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3. 6 Head to Tape Spacing Effects 

Previously, investigations have been made (see MTPR #2 and 3) 
of the effects of varying the head-to-tape spacing on the amplitude (and 20% 
pulse width) of isolated transitions. Here we report similar data applied 
to the case of repetitive transitions. The six-bit word chosen is "all ones" 

(NRZ) i. e., the record current changes sign after each bit interval. It will 
be realized that, even though an isolated pulse is insensitive to head-to- 
tape spacing because it is not rich in harmonics of wavelength shorter 
than the coating thickness, this simple conclusion need not hold in the 
repetitive pulse case. Here all the low [ <2/BPl] harmonics are forbidden 
and only the high harmonics, which are very sensitive to head-to-tape 
spacing remain. Typical data is shown in Fig. 9 , in which figure is shown 
the output voltage, at 8000 BPI, versus deep-gap field, for four cases; 

a) record and reproduce spacing = 100 p" 

b) record spacing = 100 |i", reproduce spacing = 20 p." 

c) record spacing = 20 p.", reproduce spacing = 100 p." 

d) record and reproduce spacing = 20 p." 

The curves are labeled (a^ a 2 ) where a^^ is the record-head spacing and a^ 
is the reproduce-head spacing. 

Several comments may be made about these results. 

a) Any spacing greater than the standard 20 p." at either head causes 
losses which cannot be recovered completely by adjustments of H^. 

b) Whenever the record head-to-tape spacing is large, some recovery 
of signal is possible by increasing the record current (i. e. , H^). Note that 

the (100, 20) and (100, 100) curves both peak at ~ 3500 oe, whereas (20, 20) 
and (20,100) both peak at H^ ~ 1500 oe. This recovery of signal is 
approximately 6 db, where as the loss of signal due to 100 p." record-head 
spacing is approximately 10 db. 
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Head-Tape Spacing as noted. 8000 BPI 111111 NRZ 
d = 200 pin. *^l/**2 ^^0/500 oe 

^record = 

^reprod = ^5 pin. 


Fig. 9. Output Voltage Versus at Various Head Tape Spacings 
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c) Whenever the reproduce head-to-tape spacing is large, 
a very large loss in signal occurs. This signal loss is approximately 
20 db regardless of the record current level. This is not unexpected, since 
at 8000 BPI the output signal is almost purely sinusoidal, of wavelength 
250p", and consequently the increase in reproduce spacing loss is: 

55 [(a,), - (a,),] ^ 55 [ 100 - 20] ^ ^ 

X 250 


We conclude, therefore, that at high bit densities the "uncompensated" 
record-head spacing loss is about one third the reproduce spacing loss. 

If the record current is readjusted, the "compensated" record-head spacing 
loss need only be about one fifth the reproduce spacing loss. The reproduce 
spacing loss is closely approximated by the usual 55 a/X type expression. 

3. 7 Pulse Crowding Graphs 

As called for in the contract, sufficient computer runs have been 
made to derive graphs of output voltage versus bit density. Each of the 
graphs shows the computer peak output voltage (i. e., not adjusted to yield 
absolute volts) versus deep gap field for five bit densities, those corres¬ 
ponding to 500, 1000, 2000, 5000 and 10, 000 BPI respectively. The data 
is shown in Figs. lO through 13 each figure applying to a different pair of 
heads. In Fig. 10 , the record gap is 250 (i", the reproduce gap zero, Fig. 11 
250, 75(1", Fig 12, 150, 0, and Fig. 13, 150, 75(i". The head-to-tape spacing 

on all runs was equal to 20 (i" (the nominal contact value), the tape coating 
thickness was 200 (i", the range of switching fields 100/500 oe and permeabilities 
4/2. An additional piece of information, included in all the graphs, is whether 
or not super-position (of output pulses) applies or not. Super-position only 
applies when the original recorded magnetization patterns are completely 
separate and distinct, and thus pulse crowding (or bit interaction) only occurs 
during the (linearly super-imposable) demagnetization, remagnetization, flux 

collecting and differentiating processes. 
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Fig. 10 Pulse Data Crowding 
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record 

g . = 75 uin. 
reprod 

A Superposition valid 
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Fig. 11 Pulse Data Crowding 
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Fig. 13 Pulse Cro^ 
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The following comments may be made. 

a. In all cases the effects of bit crowding do not become severe 
until above 2000 BPI. This "break point" is, with standard 
yFe 202 tape, principally dependent upon the coating thickness 
used. If standard digital 400 p." thick tape were to be sub¬ 
stituted for the present 200 p" tape, the effects of crowding 
would then become severe above 1000 BPI. 

b. Very little difference is to be noted at deep gap fields near the 
saturation level between the record gap equals 250 p" and 
150 p" data. This is, of course not unexpected. 

At low bit densities, the output voltage is governed almost 
entirely by the demagnetization effects discussed in MTPR #1. 

The details of the recording process are essentially forgotten. 

At high bit densities the recording takes place in a narrow zone 
downstream of the record gap center line (trailing edge effect). 

The geometry of this zone is not simply related to the reccrd gap 
length. It is a function of the deep gap field, the depth into the 
tape, the tape particle switching fields, as well as the gap 
length. 

The only region in which appreciable differences occur between 

the two record gap lenghts is at very high record levels (greater 

than two times saturation) and very high BPI. It will be noticed that, 

at H = 6000 oe and at 10, 000 BPI, the g . = 250 p" data is 
o ’ record 

about a factor of three (10 db) lower than in the 150 p" case. 

We may conclude that smaller gap lengths are better only in 
poorly designed regions of operation. 

c. As was mentioned in MTRP #4, at the higher BPI's, it pays 
increasing dividends to operate below the saturation (i. e., 
self-erasing) level. 
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d. Super-position is valid for all the 500 and 1000 BPI data. 

At 2000 BPI the 250 p" record gap data is not super-imposable 
above 1500 oe deep gap field. The 150 p" record gap data 
is super-imposable at 3000 oe deep gap field. 

Above 2000 BPI neither head yields super-imposable results at any 
deep gap field. 

e. The effect of reproduce gap is seen to be very slight in the 
region investigated. The zero length gap obviously will show no 
interference effects, the 75 p." case will have a null at 

X = 75 p", corresponding to 26, 666 BPI and will show negligible 
gap losses below X = 150 p" (i. e. , 13, 333 BPD. This is 
borne out in the data where even at 10, 000 BPI there exists less 
than 2 db difference. 

A practical rule emerges concerning the selection of reproduce gap 
lengths. The gap length should be as large as possible but no larger than 
the smaller of: 

a. the thickness (low BPI limit) 

b. 4/BPI (high BPI limit) 
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3. 8 Peak Shift Calculation 

In accordance with the contract, calculations of peak shift have 
been made. The computer recorded a single pair of transitions under the 
following conditions: 

d = 200 p." 
a = 20 p" 

= 100/500 oe 
= 4/2 
as noted 

as noted 

0 0 1 1 0 0 NRZ 



Percentage peak shift equals 100 ^/BPI 


Coating thickness 
Head-to-tape spacing 
Switching fields 
Permeabilities M-j/h -2 
Record gap 'j 
Reproduce gap J 
Deep Gap field 
Bit pattern 
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Table 3. 1 Percentage Peak Shift as a Function of BPI, 
Deep Gap Field, Record and Reproduce Gaps 




g . = 150 

record 

150 

250 

250 

BPI 

H 






o 

g , = 75 

reproduce 

0 

75 

0 



0. 15 

— 

— 

0. 25 

2, 000 

1500 \ 

0. 5 

0. 65 

0. 75 

0. 75 


3000 J 

1 

0. 55 

1. 6 

1. 75 


500'^ 

1. 4 

0. 9 

4. 0 

2.7 

5, 000 

1500 > 

10. 0 

9. 0 

12. 2 

11. 0 


3000 ^ 

14. 0 

13. 5 

18. 8 

17. 5 


500*^ 

12. 6 

7. 5 

23. 2 

17. 1 

10, 000 

1500 ) 

37. 7 

32. 5 

39. 1 1 

34. 2 


3000-^ 

43. 8 

39. 0 

142. 0 i 

151. 0 


Several, not unexpected, trends are evident: 

a. The percentage peak shift increases with increasing BPI. 

b. The percentage peak shift increases with increasing deep 
gap field. 

c. The percentage peak shift increases with increasing record 
gap field. 

d. The percentage peak shift increases with increasing reproduce 
gap length. 

The practical rules which emerge are: 

a. Peak shift is negligible below 2000 BPI. 

b. In the neighborhood of 5000 BPI the use of saturation* record 
current leads to peak shifts of order 10% - 20%, the precise 
value being dependent upon the gap lengths used. The peak 
shift can be rendered negligible ( < 5%) by adjusting the 
record current appropriately. 
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c. In the neighborhood of 10, 000 BPI the use of saturation* 
record currents leads to large peak shifts (>30%). The use 
of two times saturation currents causes over 100% peak 
shift. The situation may be largely corrected, however, 
by the combination of smaller record gaps (<150^.”) and 
lower record currents (about one-third of saturation). 


*Saturation means self erasing (H^ = 1500 oe for record gap = 150 p." 

(H = 2350 oe for record gap = 250 u") 
o 
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4. 0 DETECTION OF SIGNALS OF KNOWN SHAPE IN THE 
PRESENCE OF WHITE GAUSSIAN NOISE 


4. 1 Signal Representation 

4. 1.1 Restrictions 

We consider only signals 
bandwidth W. 



of finite duration T and finite 



In this case the signal s(t) is completely described by a finite set of points. 

4. 1. 1 Sampling Theorem 

Fourier integral: s(t) = 

The signal is bandlimited; the spectrum can be developed in Fourier series: 

Complex Fourier coefficients (C * : conjugate of C ): 

n n 


S(f)e>2’''' df 


c* 

n 


2W 



+j 27 !nTf 


if; 


T = 


2W 


Comparing with the Fourier integral of s(t): 

s(nT)= \ S(f)e^^’^^'^^df = 2WC* 
V J n 
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The spectrum is determined by the values s(n t), and from the spectrum the 
time-function s(t) can be computed. We need only know sampled values 

T= seconds apart to reconstruct the complete signal. Since the signal 

is of duration T, s(t) is given by a set of T/t= 2WT points . The signal could 
be represented by a point in a 2WT-dimensional space. 

4,1.3 Ideal Lowpass 

We determine the impulse response of an ideal lowpass: 



4. 1. 4 Representation of s(t) 


s(t) = S s(nT) h (t - nr) 


27 rW(t - nr) 


1 ^ 

2W 


Total number of points: 2 Wt 


2WT 
s(t) = Z 


n =1 


sin TT (2Wt-n) 
IT (2Wt-n) 


(4. 1) 


That this represents the signal s(t) is readily seen: 

-The constructed signal consitsts of a superposition of impulse responses from 

a lowpass filter of bandwidth W; hence s(t) has a bandwidth W. 

-For t = m • T = m • T,.r equation (1) becomes an identity. 

2W 
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The functions q = constitute a set of orthogonal functions; 

^ _ irx'x^TA. _\ 


ff (2Wt-n) 


< 4*00 

r r 0 n / m 

\ “ 1 i 

^ 2W 


m 


Proof: 


Parceval's equation: \ q (t) q (t) dt = \q (f) • Q * (f) dp 

v/ n in V n not 


(4.2) 


*) conjugate value 


q (t) = h{t-nT) Q (f) = H(f)* e 
n n 


j2iTfnT 


H(f): ideal lowpass 
(f) • Q * (f) df 




m 


W 


j2fff(n-m)T 
e df = 


2W 


sin ff (n - m) 
IT (n - m) 


This is for n = m, 0 else. 


4. 1. 5 Integration of Signals 


s^(t); Sj^(t) are two signals with same restrictions as stated in 1.1. 


4*00 4*00 

S s^(t) Sj(,) dt = 5 ( £ 


* <5n) ( 5s, „ dt 


n on ■“ m Im m 


(see (4. 1)) 


applying the orthogonality (2): 
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, 2WT 

S^(t)Sj(t)dL = 2W 2 V ■ "1" 

' n = I 


(4.3) 


~i. fs *sXifs,s, are vectors 
2 ^ L_o _ii' _i 

in the 2WT dimensional space: 

®o"{®ol'®02' •••®o2Wt1;' s^={s^^, s^2' •••®12Wt\ 


4. 2 Signal Detection 

4. 2. 1 Distortion of a Signal by White Gaussian Noise n(t) 

2 

n 

1 2o ^ 

_ . . . - e ; spectral density Q (f) = 

In the system signal and noise are assumed to have the same bandwidth 

— 2 " 

W: mean noise power n = W. 

*2 2 
n equals the variance o . 

Since we have white noise, n is statistically independent at different points 
in timeo The probability of a certain deformation x(t) of s(t) is easily 
computed: 



If x(t) is represented by the vector ^ = [x^' ^2' * * ’^2 WtV probability 

of the signal being in the interval 2S, 2 £ + dx is 
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p(Xj^, ^2' ”* ^2WT^ * ^1' *^2**’ ^^2WT 


p(x^, x^, ... x^^) 


I ^ > 2Wr/2 

\ 2i7N W/ 


2WT 


. e 


^ F=i 


The joint probability is the product of the individual probabilities. This holds 
for white gaussian noise. 

4.2.2 Likelihood Test 

Given an observation x(t) we have to decide which signal S j^ (t) of a given 
alphabet is the most likely one . Emitting a signal Sj^(t), the probability of 
receiving x(t) is proportional to 


p(Xi, x^, ... 1 Sj^(tj) 


= (-L- 

\27rN 


2WT/2 


W 


1 ^ , .2 

2N W ^ (4. 4) 


e o 


Assuming that each Sj^(t) is equally likely emitted, we have to choose that Sj^, 
which gives the maximum p (xj Sj^). 

Taking the logarithm of (4): 

1 , v2 1^2 2_ 

C - i.r ?(x, “ S ) - C - o»T -i.r ’••• 2/ X,S 


2N W T i ki' 
o 


2N W t ki 2N W i ki 
o o 


C is a constant, independent of index k. 

Applying equation (3); 

[c - -^T ^ ^ ®k^‘^ "^^^max =4' 
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Sj^ (t) dt = Energy Ej^ of alphabet signal 


^x(t) Sj^(t) dt = Peak output of a correlator or netched filter. 
= ^ 5 c(t) Sj^(t + t) dt; Peak = 


Decision criterion: 


max|^ ^ x(t) Sj^(t) dt - ^ 1^=^ Sj^ ( 4 . 5 ) 

4.2.3 Error Probability: k = 2 

We have to choose between two signals s (t), s,(t). To find the maximum 

o 1 

(4. 5) we may consider the difference of the two expressions 

i ^ '“t - iT • i \ ■ if 

o ^ o o 

ii = ^ x(t) (Sj-s^) dt - - (E^ - E^) 

o o 



Decision: u >0 : 

u <0 : s 
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x(t) can be 




1 ) = s^(t) + n(t) 

2 ) x^ = s^(t) + n(t) 


We make a false decision if: 

1 ) x^andu <0 

2) X and u >0 

o 


Both cases are completely symmetric; it suffices to study one, (The chance of 

having one of the two cases is 50%, because s , s, has been assumed 

o 1 

to be equally likely). 


We search for the probability of having x^^ and u < 0 


“' I S 


+ n)(s - s)dt- - (E -E), 

o 


^^®li ®li ■ ®oi^ ■ N ^^1 ■ V {(4.3)) 

N W ‘ o 


with 


As = s, - s 
1 o 

AE = E, - E 
1 o 


N W 
o 


F (=11 ^ -i) A=i - 


We know the probability density function of the n.;p (n.); we have to find 

1 n 1 

P^(u). 

u is a linear combinaUon of the random variable n^; like n^, u will therefore 
have gaussian distribution. Mean and variance for u have to be calculated: 
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__ 1 _ AP _ 

“ ' f * "i • - K~ ■•• "i = ® 

o ^ o 


u = 


^ s,. A s. - 


AE 


N W i li “ i N 
o o 


L(t) As(t) dt - ^ 
o o 


_ 2 f 1 2 2 _ _ _ 

^ ^ , because = n. * n. 


= 0 


n = N • W 
i o 


2 1 “ 2 \ 2 
o = L AS, = -zr~ I ^ s (t) dt and 


u N W 
o 


N 


J 


p(u) = 


\p% 


-,^2 (u -u)^ 
e u 


The error probability is: 

0 0 
1 


■i 


ER - ) p(u) du = 


r - ^ (u - u) 
e u 


du 


Defining erf:(iS) = 


-i-T 

2ir J 


e dx 


we find 


ER = erfcO) ; P- - ^ 

on 


N \ 

Ct 'J 


A s (t) dt 


1 

a = 
u 


A s (t) dt 


(4. 6) 
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i 2 

- 2 ^ dx (4. 8) 

I I andjAe^l can be computed from the correlation functions: 

"o“> ^0<“> = ^ "l*'* =o<‘> 

' \ V'* ^l‘‘> = 5 ’‘l“’ ®l‘‘’ 

lAeJ= l^^j(O) -^j(,(0)l; I Ae^|= l^^^(O) -fr|jj(0)l 

x^(t), Xj^(t) are the actually received signals for a "O" and a ”1" 
respectively. 

4. 3 Concerning Dropouts: Conclusions and Plan of Action 

(a) The differences in drop-out behavior are sufficiently great 
between different reels of tape to preclude any meaningful statistical 
analysis. This comment applies even if one restricts testing to samples 
of one manufacturer, or indeed, to samples slit from the same original 
doff. 

(b) The drop-out count on any one reel of tape is extremely 
dependent upon the method of testing. Track width, head contours, tape 
tension, tape speed are all factors which can alter the drop-out count 
by several orders of magnitude. 

(c) There exists evidence that, despite the foregoing, the actual 
drop-out behavior observed on some particular sample of tape being tested 
in some particular manner, is almost 100% repeatable between consecutive 
runs. By far the majority (>90%) of drop-outs occur at specific tape 
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locations and can, upon microscopic examination, be associated with 
physical defects. The number and position of such defects change only 
slowly with repeated passes of the tape. It is, therefore, reasonable 
to assume, in the computer-simulation of drop-out behavior, that the 
head-to-tape spacing is perturbed equally during both the record and 
reproduce cycles. 

(d) In fact, there is no physical distinction between drop-outs 
and the so-called surface noise. The complete noise spectrum of tape 
shows broad-band, almost white noise to which is added, at low 
frequencies, surface noise. A drop-out is nothing more or less than an 
instance of the surface noise exceeding some arbitrary level and/or 
duration. 

For the purpose of computation, however, it is a great simpli¬ 
fication to assume that the drop-outs are long in duration compared to 
the six bit word. This means that the head to tape spacing may be 
considered to be constant during the whole six bit word calculation. 
This simplification which will be adopted for the remainder of this 
programe is shown below 
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(e) Similarly, while there is no a priori physical justification 
for assuming that a drop-out extends uniformly across the whole track 
width, such will be assumed in the computer program. This is, of course, 
a considerably more restrictive assumption than that in (d) above. 
Intuitively one expects the longitudinal and transverse dimensions of a 
drop-out to be comparable. What we are really presuming here is that 

an effective head—to-tape spacing exists, in which the tape and head 
planes are exactly parallel for dimensions greater than the track width. 

The validity of this approximation will shortly be tested (see below). 

(f) We seek a typical probability density versus effective 
head to tape spacing curve. 



da 

In such a graph the probability that the head-to-tape spacing 
is between value a and a + da is p da (i. e. the shaded area). 
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Purely in the interests of simplicity, the function p(a) will 
be assumed to be an exponential of the type: 

f(a) oc e ^ ^ 

The value of the constant a. is determined in section 4. 3 
of this report. 

(g) Finally, the effect of drop-outs upon the eventual error 
rate will be calculated according to: 

ER = j ER (a) p (a) da 
20 n" 

The determination of ER (a) has been considered in the 
previous section. 

4. 4 The Head-to-Tape Spacing Probability Density Distribution 

Function 

As discussed in the previous section, we shall assume a simple 
exponential. It had been hoped that the exponent could be derived un¬ 
equivocally from existing data. Unfortunately, it transpires that only 
an upper bound can be found. 

The raw data falls into two classes, both of which were taken 
on an Ampex FR-1400 (longitudinal instrumentation recorder) at a tape 
speed of 60 ips and at a wavelength of 300 ^iin., using standard high- 
quality instrumentation tape. The first class of data, shown below in 
Table No. 4. 1 gives the number of dropouts per 100 ft of tape lasting 
more than 20 ^ec versus their amplitude. 
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The raw data falls into two classes, both of which were taken 
on an Ampex FR-1400 (longitudinal instrumentation recorder) at a tape 
speed of 60 ips and at a wavelength of 300 /xin., using standard 
high-quality instrumentation tape. The first class of data shown below 
in Table No.4.1 gives the number of dropouts per 100 ft of tape 
lasting more than 20 fisec versus their amplitude. 


Table 4. 1 


Loss of Signal 

No. of Events/lOO ft. 


CD 

p(a)da a' 

a' 

3-6 db 

181 


181 X 10“^ 40 uin. 

6-10 db 

4.8 

— 

s 4. 8 X 10 60 fiin, 

> 10 db 

0.4 

s 0. 4 X 10~ 80 ^in, 


We seek to fit this data to an exponential which is generally: 

, . -«(a - 20) 

p(a) = ae 


If we normalize so that J p (a) da = 1, then a = 8 

20 

The data in Table No. 4. 1 yields, for the three cases: 


a ^ 0. 43 
a ^ 0. 30 
a ^0.25 


The fit to an exponential is not particularly good but it is 
not felt to be worthwhile introducing more sophisticated models. Note 
also, since the original data applied to dropouts of duration greater 
than or equal to 20 fisec, we are left only with an upper bound on a. 

The second class of raw data gives the number of dropouts 
of magnitude greater than 6 db versus their duration. 
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Table 4. 2 


Duration 

No. of Events/100 ft 

.*. p(a)da 

a’ 

5-20 fjiS 

57 

710 

60 fxin. 

20-40 fiS 

6 

180 /xs 

60 liin. 

s 40 fis 

12 

- 

60 /xin. 


Whereas in the previous case the amplitude was well defined 
and the duration more or less arbitrary, here we have the converse. 
The values of a derived from Table No. 4. 2 are: 


a s 0. 25 
0 ! s 0. 33 


We conclude that the function p(a) is bounded by: 

/ \ n oc -0. 25(a-20) 
p(a) = 0 . 25 e 


The actual value of the constant a is therefore, a matter of 
some conjecture. In the hand calculation demonstrated in the following 
sections of this report a value of 0. 05 is used. The machine calculations 
of the overall error rate (section 4. 7) have used a value of 0. 08, 


4. 5 Worst Case Error Rates 

We present below tables of the error rate anticipated for seven 
levels of clipping. All the errors are assumed to be due to “noise-free 
dropouts" as defined below*. The errors are worst case in the sense 


*If the head-to-tape spacing momentarily increases sufficiently that the 
“one" signal amplitude falls below the clipping level, errors can occur 
even in the absence of noise. It is clear therefore that the choice 
of clipping level is critical and should be related to the noise level 
anticipated. 
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that the amplitude attenuation with spacing is assumed to be as large 


as is possible (i. e., the all-ones case as shown in the pulse crowding 

data graphs shown in section 

3. 7. Thus the error rates apply to a 

recording simulation using the JPL parameters. The head-to-tape 

probability density distribution 

function assumed for this illustration 

is 


pU) = 0. 05 e'"- 

Worst Case I 

Clipping Level 50% 

BPI 

Error Rate 

1000 


3000 

10-“-« 

8000 


Worst Case II 

Clipping Level 20% 

BPI 

Error Rate 

1000 

10-'-^ 

3000 

10-2-2 

8000 

10-*-“ 

Worst Case III 

Clipping Level 15% (IBM Compatible 

BPI 

Error Rate Standard Practice) 


1000 

lo "®-2 

3000 

10 -2-« 

8000 

• 

1 

o 

Worst Case IV 

Clipping Level 10% 

BPI 

Error Rate 

1000 

lO"^* ^ 

3000 

lO"^* ^ 

8000 

h-' 

o 

1 

• 

CD 
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Worst Case V Clipping at RMS Noise Level 


BPI 

1000 

3000 

8000 

Worst Case VI 

Bandpass equal to twice BPI 

Noise power density = 10“^^ 
units Vcycle 

Tape speed 10 ips (See next section) 
Error Rate 
< 10-20 

10-0 

10-0 

Clipping at RMS Noise Level 


Same as V but 


N = 10~^ unit^/cycle 


(see next section) 

BPI 

Error Rate 

1000 

10-00 

3000 

10 - 0-0 

8000 

IQ-O -0 


Worst Case VII Clipping at RMS Noise Level 

Same as V but 
N = 10~ unit /cycle 
(see next section) 


BPI 

Error Rate 

1000 

00 

o 
1—1 

3000 

10 - 2-2 

8000 

1 —• 
o 

1 

o 


Note that although worst cases V, VI and VII are related to 
the noise level, the errors computed are "noise-free dropouts" and not 
"noise dropouts". 
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It may be concluded that a simple amplitude detector employing 
a nonzero clipping level is not particularly useful at 8000 BPI. In 
section 4. 6, written by H. Lienhard, the theoretical advantages of the 
zero clippling level usable with the correlation detector are extolled. 


o The White Noise Power Spectral Density Function - 

It is necessary to establish appropriate noise levels for the 
computer error calculation. It has been decided to use three levels, re¬ 
presenting normal tape-limited background noise, one hundred times 
and one thousand times this level respectively. 


Normal tape background noise yields (cf. standard audio 
practice) a wideband (25 Kcs) tape-limited signal-to-noise ratio of about 
60 db, at a tape speed of about 10 IPS. This ratio is defined as: 


where 



wideband 


20 


€ 

VfN df 

■’ o 


€ is the signal output voltage 
f is the frequency. 


In this analysis the noise N^ is presumed to be frequency 
independent (so-called white noise) and the integral may be dispensed 
with. 


S 

N 


= 10 log 


tob 


10 


No ^ 


Now the output signal level given by the computer programs 
is about one unit (this unit and its relationship to absolute voltage is 
discussed in section 3.1. 
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Thus the value of N corresponding to tape background noise 

-11 2 ° 
is about 10 units /cycle. 

The three noise levels to be used in the error rate calculation 
-11 -9 -7 2 

are then 10 ,10 and 10 units /cycle and they correspond to wide¬ 

band (audio) signal-to-noise ratios of 60 db, 40 db and 20 db respectively. 
These three levels are chosen so that the whole range from the optimum 
down to an extreme case of amplifier noise (which pertains to low tape 
playback speeds) is simulated. 

4. 7 Signal Detection Theory 

The following questions concerning signal detection have been 
investigated: 

Influence of head-to-tape spacing (fading) 

Clocking 

Equalization of signal energies 
Comparison of amplitude and correlation detection 
Performance of correlation detection with zero decision 
level for different signal-to-noise ratios. 

All the mathematical derivations are found in the Appendix. We 
summarize the essential results: 

(a) Fading and Crowding 

The test signals for the correlation detector s (t) and s,(t) 

o 1 

are determined assuming nominal conditions (s^: "0-signal, " s^: "l-signal"): 
no fading (nominal head-to-tape spacing a^), patterns of all "ones" or 
all "zeros". The actual output signals (without additive noise) s^* (t), 

Sj^*(t) can deviate considerably from s^ and s^^ due to fading, and 
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crowding effects at high BPI. Equation (1) of the Appendix gives the 
error-probability for a false detection of the degraded signals s^*, s^^* 
in the presence of white gaussian noise. 

(b) Self-Clocking 

A big difference in the signal energies of s^(t) and Sj(t) 

means significant amplitude differences of the correlator outputs Jxs^^dt 

and fxs dt. For self-clocking detection the detectability curve should 
o 

show distinct maxima for both the "1" and the "0". With a large energy 
difference it can easily occur that, e. g. , the peaks from the s^^-corre¬ 
lator are wiped out by deteriorations of the s correlator output. Hence 

o 

for selfclocking the signal energies should be made approximately the 
same. In section 3. 0 of the Appendix it is shown that, on the average, 
the correct clocking has to coincide with a detectability peak. 

(c) Zero Decision Level 

There is another strong reason to prefer equal signal 
energies, in section 1. 0 of the Appendix we find that the decision level 
for the "1" and "0" is essentially determined by the difference in signal 
energies of s (t) and s,(t). If this decision level a is nonzero, then the 
detectability peak can fall below this level due to a large head-to- 
tape spacing and produce a false decision even in the absence of any 
additive gaussian noise. For high signal-to-noise ratios the dominant 
error source will be the head-to-tape spacing. With equal energies 
for "1" and "0" the decision level is shifted to zero (a = 0). In this 
case fading alone will not cause any errors at all. 

(d) Equalization of Signal Energies 

The highly desirable case of equal signal energies can hardly 
be realized in the tape recorder. However, we show in section 2. 0 of 
the Appendix that this case is easily approximated by attenuating one 
correlator output. Errors may be computed with formula (1) of the 
Appendix. 
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(e) Comparison of Amplitude and Correlation Detection 

By an "amplitude detector" we mean a decision scheme that 
operates on a fixed clipping level 1. The correct time for interrogation 
is supposed to be known. A "one" is detected if the signal is higher 
than I, The "peak detector", widely used in digital recording, operates 
similarly. The exact position of the peak is irrelevant for our analysis. 
Equation (3) in the Appendix gives the error probability for amplitude 
detection. The performances of the two detection principles are compared 
in section 4. 5 of the Appendix, assuming high signal-to-noise ratio 
(60 db) and low BPI (no crowding effects). With a nonzero decision 
level a, the correlation detector can be replaced by an amplitude detector 
with corresponding clipping level and signal-to-noise ratio. The decision 
errors are mainly caused by variations in head-to-tape spacing ( "drop¬ 
outs"). Since both schemes analyzed in the Appendix work on the same 

effective clipping level (« 16%), we find much the same error rates for 

“6 

both schemes (« 2 x 10 ). Digital tape machines are usually set on 

a 15% clipping level for an effective signal-to-noise ratio of approximately 
26 db. This corresponds to a 5% noise level. The major part of the 
noise is "machine noise" (crosstalk, firing from the drive motors, etc.). 

The above comparison is not valid for high BPI. If the 
neighboring bits are interacting, the output signals can suffer considerable 
d-c shifts, making amplitude detection rather impossible. (See illustration 
on following page.) 
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D-c Shift Due to Crowding 

Since the correlator detector is essentially a waveform de¬ 
tector, it can still operate properly at high BPI. 

(f) Performance of the Equalized Correlation Detector 

The full potential of the matched filter approach becomes 
apparent when the signal energies are equalized, thus yielding a zero 
decision level. Approximate computations give the following error 
rates (see section 4. 5. 1 of the Appendix), for 1000 BPI. 

Signal-to-Noise Ratio Error Probability 


60 

db 

10 -^^ 

40 

db 

10 -'« 

30 

db 

10 -'^ 

20 

db 

10-® 

10 

db 

10-5 

0 

db 

10-1 
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For 60 db signal-to-noise, the error probability is lowered by a factor 
18 

of 10 compared to amplitude detection. 

If the assumed drop-out model (section 4.2 of the Appendix) is 
essentially correct, then it has been proved that the performance of 
digital recorders can be tremendously increased by appropriate detection. 

4. 8 Computed Overall Error Rates 

This section shows the overall error rates computed by JPL4 

(see section 2.1) and JPL6 (see section 2. 2). As has been agreed 

previously, we show three graphs of error rate versus BPI for differing 

values of the steep gap fields as illustrated in Figs. 14 through 16. 

Each graph has three plots appropriate to the three noise levels (N ) 

o 

discussed in section 4. 5. The head-to-tape spacing function used is: 


200 /iin. 

150 ^in. 

75 ^in. 

100/500 oe 
4/2 

000111, frequency doubling 

The calculations pertain to the optimum type correlation detector, 
with equalization of the signal energies, as is described in sections 
1.0 and 2.0 of the Appendix of this report. 


p(a) = 0. 08 e 
Other parameters are: 


-0. 08(a-20) 


Tape coating thickness 
Record gap length 
Reproduce gap length 
Switching fields 
Permeabilities 
Bit pattern 
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Certain observations may be made concerning these graphs: 

(a) The choice of = 800 oe (see Fig. 14), that is one 
third of saturation is obviously bad. The record level is too low to 
measure against even small increases in the head-to-tape spacing. 

At a head-to-tape spacing of about 120 ^in., the 300 oe longitudinal 
field contour of the record head is just grazing the near side of the 
tape. Thus for all spacings greater than 120 ^in. ^ signal is, 
according to the modified program (see section 6. 0) written at all. 
This fact would require further thought in any continuation of this 
study. 


(b) For the case of saturation recording (H^ = 2400 oe), 
(see Fig. 15), no particular comments need be made except to 
remark that the standard IBM compatible machines have effective 
signal-to-noise ratios of about 26 db which falls about mid-way 
between the tape noise plus 20 db and plus 40 db levels. The 
overall error rates applicable to the IBM machine are expected to 
be about lO”® at 500 BPI, lO"® at 1000 BPI, lO"'^ at 2000 BPI 
and etc., while figures sound to be very reasonable. 

(c) At three times saturation level (H^ = 7200 oe), (see 
Fig. 16) we see two points of great interest. First at low BPI 
(< 1000) the error rate is decreased relative to the saturation 
level case. This is because of the ability of the head folt 

to continue to write satisfactorily on the tape even during 
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moments of long head-to-tape spacing. Some insurance 
against "drop-outs” thus occurs as is, indeed, the popularly held 
belief. Conversely at high BPI (> 1000), the high record current causes 
so much (non-linear) pulse crowding that the computer program was 
unable to detect the output correctly. Specifically the trouble appeared 
to be that "peak shift" (see section 3.8) so distorted the "detectability" 
curve (i. e.^ the difference between the moduli of the cross correlators) 
that the peak seeking algorithm used fails. This difficulty was not 
unexpected, for it has long been apparent that the use of high record 
currents at high BPI is an extremely bad situation which should be avoided. 

In the following section, proposed topics for a continuation 
of the detector study are outlined. 

4. 9 Proposed Topics for a Continued Study of the Signal 

Processing in Digital Recording 

We summarize a few topics that we feel essential for the 
optimization of the digital recording process. The main points include: 

Synchronization 

Signal design 

Noise 

Optimization with respect to dropouts 
1. Reliability in Synchronization 

a) Actual flutter data have to be incorporated in the 
calculations. 

b) The gaussian noise will randomly shift the detectability 
peaks. 

c) Crowding effects result in peak shifts. 

All three effects determine the probability of false synchronization. 
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Fig. 14 Computed Error Rate Versus BPI for Deep Gap Field = 800 oe 
(1/3 X SAT) and Three Noise Levels 
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Fig. 15 Computed Error Rate Versus BPI for Deep Gap Gield H = 2400 oe 
(1 X SAT) and Three Noise Levels 
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Fig. 16 Computed Error Rate Versus BPI for Geep Gap Field 

H = 7200 oe (3 x SAT) and Three Noise Levels 
o 
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2. Peak Seeking Procedures 

The results of point 1) will indicate how carefully the 
synchronization has to be carried out. The peak location problem can be 
treated as a parameter estimation problem. 

3. Signal Design 

The record signals for ''zero" and "one" should be designed 
to minimize crowding effects in the reproduce signal. 

4. Tape, head, and electronics noise 

Since, the actual noise power density will not be white at 
the recorder output it might be necessary to include reversible "whitening 
filters" to equalize the noise power density. This way we are able to 

find the correct optimum detector for non-white gaussian noise. 

5. Optimization with Respect to Dropouts 

So far the decision scheme has been strictly optiftiized 
only for gaussian noise. The fading, which is due to dropouts, has 
been taken into consideration only by forcing the decision level to zero.. 

A study is proposed to find proper encoding to minimize the effects 
of dropouts. 
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APPENDIX I 


1. THE CORRELATION-DETECTOR IN THE PRESENCE OF FADING 

The effect of the variations in head-to-tape spacing (a) will 
be treated as "fading." Since this fading is caused by an attenuation of the 
signals s^(t) and Sj^(t) the matched filter or correlation detector described in 
the appendix of the monthly report No. 3 still gives the appropriate solution. 
The error figures will obviously depend on the actual head-to-tape spacing. 

Assuming a non-para metric detection scheme, the matched 
filters have to be designed for some nominal spacing a^ (and, of course, 
some fixed BPI). 


We make the following distinctions; 


s (t), s,(t) are the optimum test signals for the "0" and 
o i 


"1" at nominal conditions (no fading). 


s *(t), s *(t) are the actual output signals in the absence 
o 1 


of additive noise. 


For the nominal case we find the following optimum 


decision scheme: 


yin 




©■ 


-/a. b/U 'tyJ 




u 


In the monthly report No. 3, page 15, we found the following decision 
criterion: 
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For the two signals s (t), s,(t) this means 

o 1 

or 

^ 7 f -Av/y ^ 

4dS ■= Si ft) — Sjt) 4' ^ T /'srl» 7 * 

^ fxAS ti4, i ^^ 

For the error computation we determine the distribution of the detectability u 
assuming white gaussian noise n(t): 

sr c> S^ech-O^t ^—■*>» O-^-cis. •f 

f)S^ [^yO-^n/t)Jc^^ ^ C tcxr r/" 

^ - f ^ss*o^6 


False decision: 
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Ot 


decision level = 





, Uj^ are the actual detectability peaks in the absence of noise 




= error rate for fixed spacing a 


P = P(a); q = q(a). 


2 . EQUALIZATION OF SIGNAL ENERGIES 

We modify our detector to obtain equal detectability peaks 
for "0" and "1" in the nominal case (no fading, no noise). 
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For equalization: 


Decision; (o^= 0) 


U - 



' 3, V/ 

* ^ 

^ ^ o 


o 

JS=p- 

c 

- 


/s.V/ 

tt 4 
/ 


# 4 

o 



0 


Distribution of u; 

We find in the same way as in the preceding paragraph 


^ Js, /or =0 


If the signalss (t) and s, (t) v/ere actually designed to give equal energies, 

® 2 \ 
It 2 t> 2 2 (t 2 , , ^2 _2 f 1 + k \ 

J s dt = J s dt, we would find o = N J s dt, hence a = o 2 J 


Error probability: 



RR 67-36 


82 





AMPEX 



u being the detectability peaks in the absense of noise. 


For equal energies s^(t) and Sj^(t) we would find 

((^1 ; 


? 




u vJ- 


T’ 


OToT S \ -(2) 


Since o comes close to a^, our attenuator scheme and the real equal energy 
scheme give much the same error rates. Thus for error-probability calculations 
formula (1) may be used. 

3. EXTRACTION OF TIMING INFORMATION 

Our selfclocking code should provide enough information to 
control the clocking. We give a short verification of the fact that the detectability 
peaks indicate, on the average, the correct time of the occurrence of either s^ or 
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u being the output of a correlator: 

U.it) - ^ K W S {x-t)e^v ■ XLX)- S CtJ ■•■rii’t) 


Averaging over t: 




Jslr) s LT~t}^r +- 
f S(v) s(r~t)^T 


n It) S (t- t) 


Jr 


jr o 


Since StT:) — -S (T - ir') It follows that 

(iteWT --2- |smccr-tWr 

or 

fs\T)^T ^ fsa)s[T-t)^r 

Hence 


^lo) UU) 


t = 0 indicates the occurrence of s(t) and we have proved that at t=0 the 
correlator output must have, on the average, a maximum. 

4. AMPLITUDE AND CORRELATION DETECTION IN THE PRESENCE 

OF FADING AND WHITE GAUSSIAN NOISE 
4.1 Amplitude Detection: 
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signal s(t) 
t = decision level 
x(t) = n(t) no signal 
or 

x(t) = s(t) + n(t)i=|?signal present 



Decision: 


xCt©) ^ t i=4> signal present 
s(to) < t 1 ==!? no signal 


This would constitute an amplitude detector for NRZ - Code. 
Noise n(t): HIf) ^ a ^ ) 


at time t 

o 


X ^ S tvU,) <3nr 


fx ' pn ( <=T- iV 


False decision: 



ay 

~ pni^) 
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^a> 


ft ^ 1 










N = RMS - noise voltage s' “ / ^ 


4,2 Fading Model 

Since the output voltage of a tape recorder decreases, in a first 
approximation, exponentially with the spacing a, we define the fading: 


Sit) » • e 






(4) 


S : voltage for nominal spacing a = a 
n o 


(f=o) 


k « 0.134 for 1000 BPI 


4.3 



Distribution of Head-to-Tape Spacing 




e 


^ /f) = ‘S 


(5) 
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The average error rate is 

^ J fs (^) - 

c> 

For the following computations we shall assume: 

c = 0.05; a = 20 M" 
o 

p (^ ) = e ^ 

4.4 Correlation -Detector 

To extimate the performance of the correlation detector, we 
must make certain assumptions about the waveforms of s^(t) and Sj^(t). 



Decision level a = |- ^Jsj^^dt - J s^dt^ = ^ T 

* — Icf k£ 

Attenuated signals: s, =s,e s* = se 

11 o o 




^ J'as S*^ 


U. r 


As 


/7 ^ ’H' ^ AX 

^ Ts 3 ^ - Ts c. 

J 
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-oC 


2 T t'- r / 


*r- r 

' ^ I 


e. — 


r —/r 


^ o{ ^ TJ'^^ ^ %7 y /V".' 


< « i 




cs. /5^ 


W[^'‘’-'Aht 


>, A , A* 

For the error-rate we found in Paragraph 1) 

^ ■ f flP') * T 


# /StT 


Comparison of Amplitude Detection and Correlation Detection 


Amplitude detection: 




A/f j 


S «• -S*« € J 




For both detection schemes* 


7i '- i. f /P/) + T flp^) 
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Estimate of time - bandwidth product: 



(sampling theorem) 


Choose 




3 

For signal -to-noise ratio we choose ■y= 10 (60 dB). If y is large, 

becomes negligible. Consequently: 





Amplitude Detection 


Conclusion: 


Correlation detection 





To analyze the correlator detector , it can be replaced by an amplitude detecto r 
operating under a different signal-to-noise ratio y This is only true for 
high y and nonzero decision level a for the correlation detection. 

Rough estimate of error rates: 



clippling level 
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R: Uncertainty region due to gaussian noise R - 


False detection if S falls below ^ -f N 


A 

S ^ -A ^ 

s g F ^ r ^ = s 




To the attenuation — corresponds a critical spacing i , every 4 ^ 

S k k 

give rise to false detection: 


Oo 

ER‘ J fk/7) oi^ ‘ e: 


(see 4.3) 


With krO.134, a =20/i", c = 0.05, we find 

o 


13.38 


ER (Amplitude Detection) = e 


— 13 32 

ER (Correlation Detector) = e 


(« 2- 10"®) 


Both schemes give almost the same error rate. This has to be the case, since 
both detectors have the same effective clipping level (approximately 16%) and 
at this clipping level the gaussian noise has a negligible effect. 


4.5.1 


Zero decision level for the correlation detection 


We assume equal signal energies for s^(t) and Sj^(t): 

2 2 

fs dt = [ s, dt . If this is not apriori the case, we come very close to 
*' o "1 

it by applying the equalization scheme of paragraph 2). Our decision level 
becomes zero (a= 0). With the assumptions of paragraph 4.4 we easily 


find: 
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/f.V/' - /r,"^ - 

4 = - Ts^ ; f 

^ - / — ' 

^ ^ ^ - livT e. 




Total error rate (over all spacings 



We approximate 


I 

- A 

by £ 




CA,=| 


/ ( 3 

defining «*iwrr’'c us'.io'’ 



6 
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We change the variable of integration: 


O ^l(z) 




Z I I ^1 


Z - / : I - ^ 


- > z^ 


We approximate l(Z) by a straight line through (0, 1) and (z , —) 


-dcZ, 

C. “ « 




/.ST - /fi*' 

” -e/T 


f-c^' k ^ ^/doo-BT^^ Cc^'^o ^ ^33 






Linearization: 


Jl(z) ^% 


*»» * I 

.. 


/ - % 


== / -r<f 


^,7f • 




r~2 U 
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Comparing this to the case of nonzero decision level (see 4.5) we see that 

18 

the error rate is lowered by almost a factor of 10 for similar signal energies 
and equal signal-to-noise ratios. 

Different signal-to-noise ratios: 
define =20 log y 

log Z =-0.373 (10 log 1.5 + fork = 0.134 

o dB 

10 log 1.5 = 1.761 


- logZ «Z ER = 0.4.Z 

dB o o o 


60 dB 

-23 

CO 

CM 

1 

o 

4 • lO"^^ 

40 dB 

-15.5 

3 . lO"^^ 

1.2 • lO"^® 

30 dB 

-11.85 

1.4 . 10"^^ 

-12 

0.56 • 10 

20 dB 

. -8.1 

1/1.25-10"® 

3.2 • 10"^ 

10 dB 

-4.4 

• 

0 

1—' 

O 

1 

1.6 • 10"^ 

1 

0 dB 


2.2 • 10'^ 

1 

o 

1—1 

• 

00 

00 

• 
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APPENDIX II 


1. INTRODUCTION AND SUMMARY 

In conformance with the extension of the contract, a short 
program of experimental validation was undertaken at the conclusion of the 
original study. The validation sought is concerned with the basic tape 
recorder simulation (JPL-4) and is in no way related to errors and reliability. 

To this end, appropriate recording apparatus was assembled and 
actual output waveforms photographed. These were then compared with the 
computed waveforms. Although the extensibn contract calls for a set of 
four such comparisons, in fact, eighteen runs were made in order that a 
greater range of prarmeters might be covered. 

In every case, the agreement is better than had been anticipated 
earlier. In the case of saturation level recording the program is found fo be 
adequate up to at least 2500 bpi, frequency doubling (i. e. 5000 bpi NRZ). 

At lower levels the program is usable without serious error up to bit densities 
of 7500 bpi (i. e. 15, 000 bpi NRZ). 

It is concluded that, since there exists no justification for 
operating at the saturation level at high bit densities, the existing program 
(JPL-4) is perfectly adequate for all foreseeable engineering purposes. 
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2. DESCRIPTION OF THE DIGITAL RECORDING APPARATUS 

A DATAPULSE output unit 903 is connected to a specially made 
bistable digital record amplifier. The output is connected directly to an 
Ampex FR-1400, 11 turn per leg, 150 microinch gap length, 3 mil lamination 
Moly-permalloy record head. The amplifier is able to drive +125 mA current 
through the turns on one leg which may be computed to result in a deep 
gap field of 4400 oersteds. For a 150 microinch gap length the saturating 
field (i. e. that which puts the 500 oersted contour at a distance of 420 
microinch above the head plane) is about 4500 oersteds. The magnitude of 
record current is thus adequate, a conclusion borne out by the ability of 
the apparatus to self-erase previously recorded information. 

The tapes used are Ampex 771 and Ampex 772, of coating 
thickness 400 and 200 microinches respectively, both being 35% by volume 
■yFe_0 instrumentation tape on 1 mil Mylar. The maximum remanence is 
about 1350 guass, the intrinsic coercive force (M = 0) is 280 oersteads, giving a 
linearized demagnetizing permeability of about 5. The remagnetizing per¬ 
meability may be assumed to be about 2. The range of switching fields is 
approximately 100 - 500 oersteds. 

The reproduce head is an Ampex FR-1600 with 35 turns per leg, 

25 microinch gap length, 3 mil lamination of Moly-permalloy which is connected 
to a field effect transistor pre-amplifier which provides +40 db gain, flat 
within 0. 5 db between 1 Kcs and 2. 0 Mcs. 

The tape transport is a precision type, with vacuum column tape 
tension control and a lightweight capstan assembly with servo control. The 
machine provides 240, 120, 60, 30, 15, 7. 5, 3. 75, 1. 875 ips tape speeds, 
forward and reverse. 


RR 67-36 


96 




AMPEX 


The record current, measured on the head terminals, shows a 
negligible rise time of less than 0. 5 microsecond. The two frequencies 
proposed for use correspond to durations of 13. 5 and 27 microseconds, viz: 



At a tape speed of 60 ips this yields 625 bpi, at 30 ips 1250 bpi, and so 
on through 3. 75 ips yielding 10, 000 bpi (i. e. 20, 000 frpi maximum). 


The 8 bit frequency doubling word used throughout the whole 
validation program is 


0 10 10 0 11 


3. COMPUTER INPUT PARAMETERS 


The parameters used in the computer simulation correspond as 
closely as possible to the experimental situation. They are: 


Head-to-tape spacing, a. 
Tape coating thickness, d. 
Record gap length 
Reproduce gap length 
Deep gap field, H^, 


20 |j.in. 

200 or 400 pin. as appropriate 


150 pin. 

25 pin. 

As appropriate (calculated on 
the basis 35 oe per mA of head 
current) 
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Tape permeabilities, = 5/2 

Tape switching fields = 100/500 oe. 

Bit densities = 625, 1250, 2500, 5000 

and 10, 000 bpi 
(as appropriate) 

In order that the experimental waveforms could be "inputed" 
to the computer, the following procedure was used: 

a) The original Polaroid oscilloscope picture was enlarged 
400% and made negative. 

b) The amplitude of the waveform was measured, by hand, 

at some 60-70 points per 8 bit word. This corresponds to a point density of 
about 40 per inch on the original photograph (and to about 8 points per bit 
length). 


It may be seen by comparison of the experimental and reconsti¬ 
tuted waveforms that this procedure is sufficiently accurate to reproduce all 
the essential features of the waveforms. 

In addition to computing the waveform and comparing it with 
the experimental waveform, the computer was also programmed to compute the 
normalized cross-correlation between the two waveforms. This function, 
considered a "figure of merit" by JPL^ is: 


e(t) c(t + T) 


dt 


F.O.M. = 


J[e(t)j^ dtj[c(t)] 
T T 


dt 


V2 
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where e (t) and c(t) are the experimental and computed waveforms, of 
respectively, both of period T, This FOM is plotted in each case for 
values of T between O and T. 

4. RESULTS 

For each test we show three figures. The first is a high resolu¬ 
tion photograph of the experimental waveform. The second is computer graph 
showing the computed and reconstituted experimental waveforms (i. e. compute 


versus measure.) The third is the cross correlation FOM graph. 

The experimental conditions are: 

Test No. Coating Thickness Deep Gap Field 

(microinch) (oersted) 

BPI 

1 

400 

4500*' 


625 

2 

400 

4500 

> 1 X sat. 

1250 

3 

400 

4500 

2500 

4 

400 

4500^ 


5000 

5 

400 

2250 


625 

6 

400 

2250 


1250 

7 

400 

2250 

}V2 X sat. 

2500 

8 

400 

2250 


5000 

9 

400 

2250 


7500 

10 

400 

2250^ 


10,000 

11 

400 

1500 ^ 
( 

1 

► 1/3 X sat. 

7500 

12 

400 

1500 > 

1 

10,000 

13 

400 

1100 1 

j 

1 1/4 X sat 

7500 

14 
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15 

200 

2500^ 



6?5 

16 

200 

2500 


1 X sat. 

1250 

17 

200 

2500 



2500 

18 

200 

2500^ 



5000 
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Fig. 1 Coating Thickness 400 jo. in. lA 

Deep Gap Field 4500 oe (1 x saturation), 625 BPI 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross-Correlation F. O. M 
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Fig. 2 


Coating Thickness 400 ixin. 

Deep Gap Field 4500 oe (1 x saturation), 1250 BPI 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross-Correlation F. O. M. 


2A 





















CT> 

'<i 

I 

(jJ 

(T> 


O 

Cn 



3A 

Fig. 3 Coating Thickness 400 [lin. 

Deep Gap Field 4500 oe (1 x Saturation), 2500 BPI 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross-Correlation F.O, M 
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Deep Gap Field 4500 oe (1 x Saturation), 5000 BPI 
a) Experimental Waveform 


b) Computed versus Experimental 


c) Computed Cross-Correlation F. O. M 
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Fig. 6 Coating Thickness 400 ixin. 

Deep Gap Field 2250 oe (1/2 x Saturation), 125 0 BPI 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross-Correlation F. O. M 
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Fig. 7 Coating Thickness 400 jiin. 

Deep Gap Field 4500 oe (1/2 x Saturation), 2500 BPI 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross—Correlation F, O. M. 
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Fig. 8 Coating Thickness 400 fiin. 

Deep Gap Field 2250 oe (1/2 x Saturation) 5000 BPI 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross-Correlation F, O. M. 
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Fig. 5 Coating Thickness 400 jxin. 5 A 

Deep Gap Field 2250 oe (1/2 x Saturation) 625 BPI 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross-Correlation F. O. M 
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Fig. 9 Coating Thickness 400 |i.in. 

Deep Gap Field 2250 oe (1/2 x Saturation), 7500 BPI 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross-Correlation F. O. M. 
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Fig, 10 Coating Thickness 400 (xin, 10 A 

Deep Gap Field 2250 oe (1/2 x Saturation) 10, 000 BPI 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross-Correlation F. O. M 
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Fig. 13 Coating Thickness 400 fxin. 13A 

Deep Gap Field 1100 oe (1/4 x Saturation) 7500 BPI 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross-Correlation F. O. M. 
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Coating Thickness 400 |ain. 

Deep Far Field 1100 oe (1/4 x Saturation) 10, 000 BPI 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross-Correlation F. O. M. 
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Fig. 15 Coating Thickness 200 |J-in. 

Deep Gap Field 2500 oe (1 x Saturation), 625 BPl 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross-Correlation F. O. M. 
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Deep Far Field 2500 oe (1 x saturation) 1250 BPI 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross-Correlation F. O, M. 
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Fig. 17 Coating Thickness 200 hiin. 

Deep Gap Field 2500 oe (1 x Saturation) 2500 BPI 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross-Correlation F. O. M. 
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Fig. 18 Coating Thickness 200 fiin. 

Deep Gap Field 2500 oe (1 x Saturation), 5000 BPI 

a) Experimental Waveform 

b) Computed versus Experimental 

c) Computed Cross-Correlation F. O, M, 
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